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■ We develop a numerical method for detecting a vortex unbinding transition in a two-dimensional 
' system by measuring large scale fluctuations in the total vortex dipole moment P of the system. 

»' ' These are characterized by a quantity T which measures the number of configurations in a simulation 

for which the either P^ or Py is half the system size. It is shown that T tends to a non-vanishing 
{3J[)' constant for large system sizes in the unbound phase, and vanishes in the bound phase. The method 

^ , is applied to the XY model both in the absence and presence of a magnetic field. In the latter case, 

■ the system size dependence of J- suggests that there exist three distinct phases, one unbound vortex 
phase, a logarithmically bound phase, and a linearly bound phase. 

CN ! PACS numbers; 05.10.-a, 64.60.-i, 75.10.Hk 
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O ' Introduction - Topological defects play a crucial role in two dimensional classical systems pt] and in 1+1 dimensional 

quantum systems The paradigm of these is the XY model, which is known to undergo a vortex- unbinding transition 
in the Kosterlitz-Thouless (KT) universality class Highly analogous transitions occur for vortices in superfluids 
^Jj \ and thin- film superconductors, as well as for dislocations and disclinations in two dimensional crystals. Wordline 
' dislocations in one dimensional quantum systems are important for describing tunneling Q, and the question of 
c/2 whether these are bound or unbound is closely related to whether the system is metallic or insulating. 

The KT phenomenology has been highly successful in describing defect unbinding in a variety of situations. For 
the XY model, fitting to the expected finite-size scaling of the helicity modulus in simulations yields an estimate 
of Tkt — 0.892 However, the concepts of vortices and defect unbinding are more general than the systems 
in which a KT transition takes place. In some situations a symmetry-breaking field, for example, a magnetic field 
tending to ahgn the spins in the XY model, may be present A recent renormahzation group (RG) study @] 
Q ■ has suggested that vortex unbinding does occur in such systems, although the transition is considerably altered from 
O ] the KT behavior. It is thus important to develop criteria by which one may determine whether a given system is in 
a bound or unbound vortex state which are independent of precise matching to the KT theory. In this article, we 
I ' present a method by which this may be accomplished, which focuses on a measure T of extreme fluctuations in the 
^ , system vortex dipole moment. Using a Langevin dynamics simulation of the XY model, we show that the method 
' locates Tkt with reasonable accuracy. We then include a magnetic field in the Hamiltonian, and show that there are 
^-H ' both bound and unbound vortex phases. The bound phase has two distinct behaviors: for smaller fields, T vanishes 
\l ' as a power of the system size L; for larger fields, T vanishes exponentially with L. The two behaviors are consistent 
' with the results of the RG analysis 0, which predicted both a logarithmically bound vortex phase and a linearly 
1 confined one, in addition to the unbound (deconfined) vortex phase. 

■ Characterizing Vortex Dipole Fluctuations - A total vortex dipole moment may be defined for configurations of any 
-j—; . system containing vortices with well-defined locations. For concreteness, we will work with the XY model (planar 

■ spins of fixed length) on a square lattice. Let l^Oij denote the angular difference between nearest neighbor spins ij, 
^ , which we reduce to the interval — tt < A^ij < tt by adding or subtracting 27r if necessary. The vorticity qi around an 

elementary plaquettc P is then qi = (l/27r) AOij, which takes on the values —1, 0, 1 For any configuration of 

the XY system this rule assigns vortex charges to the sites Ri of the dual lattice. The corresponding dipole moment 
could be defined by 
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P = Y,q^R.. (1) 



However, there is a problem associated with the use of periodic boundary conditions: for a. L x L system one may 
^ [ add {nfL,n\L) to Ri (with nf , n\ integers) and retain a perfectly sensible definition of P. For a given configuration, 
■ - ' we can use this to reduce P so that its components are restricted to the interval — L/2 < P^^y < L/2; or we can 
extend its definition by adding factors of L so that P remains a continuous function as a vortex crosses a periodicity 
boundary. We will refer to these alternate definitions of the dipole moment as Pred and Pea:t , respectively. Note that 
^red ^'^^^ jump discontinuously by L whenever a vortex crosses a boundary. 

The bound vortex phase has the property that the dipole moment remains finite, since the pairs do not separate; 
in the unbound phase, the diffusion of vortices implies a diffusion of the dipole moment, so that its magnitude can 
become arbitrarily large. In a finite system, this presumably translates into the statement that the diffusion of P^xt 
becomes very small below the unbinding temperature. 
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A feature of vortex unbinding transitions is that the transition occurs via rare but extreme fluctuations |g] . One way 
to characterize such fluctuations is to look for system configurations with extreme values of Pred- Loosely speaking, 
if we wish to characterize a configuration in terms of an effective single vortex-antivortex pair, when P^ or Py is L/2, 
the pair is at its maximum separation. If such extreme configurations persist as L — > oo, the system is then in an 
unbound phase. By contrast, we expect the number of such configurations to vanish in the large size limit when the 
system is in a bound phase. 

This expectation may be more carefully justified by considering effective theories of vortices in their bound and 
unbound phases. An unbound vortex system behaves as a two-dimensional metal in which the discreteness of the 
underlying vortex charges may be neglected. An effective Hamiltonian takes the form jl^ 

Kil = f / d'xd^x'pix)G{\x-x'\)pix') + d'xp^ix), (2) 

where G{r) ~ — (27r)~^ log(r/a)+const, with a a microscopic length of order the lattice constant of the system, and 
p is an effective chemical potential, essentially the core energy of the vortices. For this continuum system, we can 
adopt a definition of the (x-component of the) dipole moment which incorporates the periodic boundary condition, 

P^ff =i,Jd^xp{x)sin[^^] 

= ^Im{p{k, = 27r/L, ky = 0)}. (3) 

The fraction J- of configurations with P^f-f = L/2 is thus given by the probability of finding Im{p(kx = 27r/L, ky = 

0) = TT. This is easily evaluated by re-expressing H^^l in terms of a wavevector sum instead of a real space integral. 
Noting that the Fourier transform of G{r) for small wavevectors is G{k) ~ one obtains 

J. e^I^[^^(^)'+^l « (4) 

for the probability of obtaining an extreme dipole fluctuation in the large size limit. Notice the does not vanish as 
L oo, supporting our argument above that large dipole fluctuations survive in the unbound vortex state. 
To analyze the bound vortex state, we focus on the two dimensional Coulomb particle Hamiltonian 

HcG = Y J d^^d^'^'M^)G{\x ~ x'\)m{x'), (5) 

where m{x) is an integer degree of freedom, and the partition function involves a sum over all complexions of m 
satisfying J d'^xm{x) = 0. The probability of an extreme fluctuation in the system dipole moment may be expressed 
as 

oc ^ e-"'^'=^^"'^/'^d[Im m{kd) - n] 

{m} 

dXe-'^"" g-ffcG[m]/TgiA/m m(fc,)^ 

{m} 

where we have adopted the same definition of Pi^ff as in Eq.^, with p{x) replaced by m{x), and kd = (27r/L,0). In 
the bound state the integer nature of m{x) may not be ignored; however, we can make progress by adopting the dual 
description of the partition function . This involves employing the Poisson resummation formula to rewrite Eq. ^ 
as 

J-OO J 

The integration over the continuous field 4> may be carried through, with the result 

{"} 

A bound vortex fixed point is generated by replacing K with a renormalized value, and exchanging the sum over 
integers n in Eq. N with a functional integral over continuous fields [nil . The resulting model represents the rough 
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phase of a solid-on-solid model. Once the integer sum has been replaced by an integral, the \Im n{kd) term in the 
integrand may be shifted away and the functional integral in fact has no dependence on A. It then immediately follows 
that !F = ioT the bound vortex phase. 

These considerations lead us to expect that one should observe large vortex dipole fluctuations in the unbound 
phase, but not in the bound one. We now demonstrate this is indeed the case using a Langevin dynamics simulation. 

Simulation - Our simulations focus on the XY model for which we assign dynamics to the spins and coupling to a 
heat bath to generate a distribution of configurations. The equations of motion for our system are taken to be 

dt^ " se{x) ^ dt ■ 

r is an effective moment of inertia for the XY spins, which for simplicity we set to 1 in the simulations, C is a random 
torque which models coupling to a heat bath, and 77 is a viscosity. To satisfy the fluctuation-dissipation theorem, the 
random torques are drawn from a distribution satisfying < ^(r, t)Q{r' , t') >— 2rjT5p^pi5{t — t') with T the temperature 
of the system. Finally, our XY Hamiltonian is 

HxY^-K ^ cos[e{x) -eix'y^-h'^coseix) (lo) 

where we take the angles Q{x) to reside on an L x i square lattice. To perform the simulation, we have discretized 
the time derivative in Eq. ^ and used a standard random number generator to generate a realization of C,{r, t) 
at each time step. A typical run consists of 10® Langevin sweeps for equilibration, followed by 9 x 10® measurement 
steps. In accumulating the data, we repeated runs for each set of parameters with ~ 10 different seeds, allowing us 
to estimate our statistical errors. Simulations were performed for system sizes as large as L = 199, although most of 
the simulations were in the range 19 < L < 59. 

Our measurement consists of counting the number of times a component of the system dipole moment Pext passes 
through [n + for n any integer. We then plot the number of such events divided by the total simulation time, 
yielding a measure of the fraction of configurations T for which the system has attained its maximal value. One 
advantage of the Langevin dynamics approach is that the vortex dipole moment Pext changes by several steps with 
each Langevin sweep, but these steps are always much smaller than L except for very small values of L. This allows 
us to detect when Pext has passed through (n -I- i)L even if in the immediate time step before and the step after Pext 
was not measured to be precisely at this value. A larger number of events can then be accumulated than one might 
in a Monte Carlo simulation employing a cluster algorithm, since the configurations generated in the latter are not 
related in any simple way, forcing one to count only configurations for which Pext is precisely {n -\- \)L. Note that we 
count passages in both directions; this tells us how often Pred visits its extremal value. 

As a check on the method, we first present results of simulations in the absence of the symmetry-breaking field, for 
which vortices unbind in a Kosterlitz-Thouless transition. Fig. 1 illustrates these for L = 59, h — Q and K — 1. One 
can see J- decreases sharply as the temperature approaches T k, 0.9 from above, so that J- appears to be vanishing 
quite close to the accepted value of = 0.892. The differing behavior of T above and below the transition can be 
further confirmed by examining its size dependence: for T below the transition temperature, T decreases with system 
size, apparently approaching zero as L — > 00, whereas for higher values of T, T increases, approaching a constant 
value from below. These differing size dependences strongly support the idea that T distinguishes the bound and 
unbound phases. 

We now turn to the case /i > 0, for which we show some typical results in Fig. 2. The behavior of as a function 
of system size takes three differing forms, depending on the values of h and T. At high temperature and low fields, 
J- clearly approaches a non- vanishing value in the large system size limit (e.g., h — 0.18 curve in Fig. 2). Unlike 
the /i = case the asymptotic value of J- is approached from above, indicating the that symmetry breaking field has 
had some effect, although for small enough h the vortices remain unbound. At intermediate values of ft,, we find T 
decreasing as a power law in 1/i {h — 0.3 curve), indicating a bound vortex phase behaving as one would expect for 
logarithmically interacting vortices. Finally, at the largest values oi h {h — 0.5 curve), J- decreases exponentially with 
L, as might be expected for vortices interacting with a linear binding potential. 

Our results indicate that in the presence of a magnetic field, there exists an unbound vortex state and two different 
bound vortex states in the XY model. This behavior is precisely what was predicted in the RG analysis of Ref. . At 
low temperatures, vortex-antivortex pairs are connected by a string of overturned spins, leading to linear confinement. 
As temperature is increased, fluctuations can lead to a roughening of the strings binding the vortices, driving the 
effective string tension to zero. The vortex-antivortex pairs nevertheless retain a logarithmic attraction and remain 
bound as the system passes through the transition. At still higher temperatures, a second transition occurs in which 
the vortices do ultimately unbind. A remarkable feature of these transitions is that they are extremely continuous 



3 



, so that no singularity is expected in the free energy as the transitions occur. Measurements of the specific heat 
and magnetization in our simulations show no indications of any such singularities. Thus, one is forced to probe the 
vortices directly in order to detect their unbinding, as we have done by probing the system vortex dipole moment. It 
is interesting to note that, in the absence of singularities in the free-energy, it is unclear whether vortex unbinding 
should be thought of as a true thermodynamic phase transition. To our knowledge, this is the first example of a 
defect-unbinding transition that docs not have the full character of a phase transition. 

In summary, we have developed a method of detecting vortex unbinding in two-dimensional systems by tracking 
extreme fluctuations in the vortex dipole moment of the system. For a system of size linear L, the fraction of 
configurations T having a dipole moment component Px or Py equal to its largest value consistent with periodic 
boundary conditions {L/i) approaches a constant for large sizes when vortices are unbound, and vanishes when 
vortices are bound. We have demonstrated this for the Kosterlitz-Thouless transition in XY model, for which !F 
can locate the transition temperature with reasonable accuracy. In the presence of a magnetic field tending to orient 
the spins, we have found that there is an unbound vortex state and two possible bound vortex states, one consistent 
with logarithmic binding of the vortices, the other with linear confinement. The method presented here is easily 
generalizable, and should be applicable to systems with other topological, point-like defects. 
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FIG. 1. F vs. temperature for XY system with ft = 0, L = 59. F drops sharply in vicinity of the known Kosterhtz-Thouless 
temperature. Left inset: F vs. i<x T = 0.8, demonstrating F vanishes in the bound vortex state. Right inset: F vs. L 

for T = 1.2, demonstrating F approaches a non-zero constant in the unbound vortex state. 
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FIG. 2. F vs. system size for XY system for h = 0.18, 0.30, 0.50 for fixed temperature T — 1.2 F approaches a constant for 
large system sizes for smallest field, vanishes as a power law for intermediate field (straight line on log- log plot), and vanishes 
exponentially for largest field (downward curvature on log- log plot). The three different behaviors indicate one unbound vortex 
phase, logarithmically bound phase, and a linearly confined phase. 
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